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PART  I.  This  report  presents  results  of  continuing  studies  of  the  Coral  Pink  Sand  Dunes  Tiger 
Beetle  funded  by  a Challenge  Cost  Share  Agreement  between  the  Bureau  of  Land  Management 
and  Randolph-Macon  College.  Included  in  this  report  is  new  information  obtained  in  the  year 
2008  and  an  updated  population  viability  analysis.  The  total  adult  population  size  estimate 
of  C.  albissima  for  2008  was  1072,  with  a 95%  confidence  range  of  959  to  1185.  This 
count  represents  a significant  increase  from  the  700  in  2007  but  a decline  from  the 
previous  even  year  cohort  estimate  of  1112  in  2006.  The  2008  estimate  also  indicates  a 
6-year  period  of  counts  of  1112  or  less  since  the  highest  ever  count  of 2944  was  recorded 
in  2002.  This  6-year  period  also  includes  all  but  one  of  the  lowest  counts  since  counts 
were  initiated  in  1992.  The  densities  and  distribution  pattern  among  the  individual 
swales  were  generally  similar  to  that  in  most  previous  years.  Highest  adult  numbers 
were  in  swales  JK  with  280  adults,  2-3  with  1 71,  H-WH with  161,  1 with  118,  public  with 
63  and  QR/EQR  with  107.  As  has  been  the  norm  in  low  abundance  years,  numbers  of 
adults  in  swales  outside  of  the  primary  core  area  were  very  low:  only  15  adults  were 
found  in  all  swales  north  of  swale  4 and  only  11  adults  were  in  the  D swales.  Only  two 
adults  were  counted  in  the  AAA-AA  area  at  the  far  northern  end  of  the  dune  field,  but 
larval  numbers  in  these  northern  swales  were  higher  (182  total)  than  in  other  recent 
years.  The  total  percent  of  the  population  within  the  protected  swales  (within  the 
conservation  area)  was  80%,  continuing  the  pattern  seen  in  all  years  since 
implementation  of  the  conservation  agreement.  Continued  monitoring  of  rainfall  totals  at 
the  Kanab  station  has  continued  because  of  the  importance  of  rainfall  in  driving 
population  trends.  The  rainfall  amounts,  both  total  for  April  through  October  (the  period 
of  adult  and  larval  activity)  and  April-June  (the  critical  period  for  oviposition  and 
recruitment)  were  the  lowest  ever  in  2008,  2.82  inches  and  0.22,  respectively.  This 
compares  with  the  previous  low  April-October  total  of  5.15  in  1993  and  the  previous  low 
April-June  total  of  0.41  in  2002.  It  is  also  notable  that  the  totals  for  2006,  2007  and 
2008  were  significantly  lower  than  any  other  three-year  period,  and  consequently  project 
continuing  low  adult  numbers  for  at  least  the  next  two  years. 

PART  2.  The  new  PVA  is  an  updated  version  of  the  1999  one  and  considered  to 
be  more  reliable  because  it  includes  more  years  of  data  (10  versus  7)  and  because  the 
population  estimates  in  these  past  10  years  are  much  more  accurate.  The  primary  result 
of  this  PVA  is  that  the  CPSD  Tiger  Beetle  continues  to  be  at  risk  of  extinction.  If  the  past 
10  years  parameters  are  used  (srowth  rate  of 1.058,  standard  deviation  of  srowth 
rate=stochasticitv  of  .646,  carryim  capacity  of 3500  adults,  startim  population  of 1200). 

the  probability  of  the  population  becoming  extinct  is  .32  and  the  probability  of  the 
population  declining  to  50  individuals  is  .87.  However,  substituting  a growth  rate  of  1.35 
while  keeping  the  other  values  the  same  reduces  the  extinction  probability  to  .001  and  to 
50  individuals  to  .003.  This  comparison  demonstrates  the  major  effect  of  growth  rate  on 
extinction  risk.  The  growth  rate  over  the  10  years  was  1.058,  based  on  positive  growth 
rates  in  only  three  years  (1.25,  1.47,  2.32);  all  other  years  had  negative  rates,  thus 
producing  a low  mean  growth  rate.  The  only  other  parameter  which  had  a significant 
effect  on  the  resulting  extinction  risk  is  stochasticity.  If  this  value  is  set  at  .33,  extinction 
probability  is  very  low  in  most  scenarios,  but  if  it  is  increased  to  .95  (a  value  that  was 
seen  among  some  of  the  years)  extinction  risk  increases  to  near  1.0.  These  results  suggest 
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that  long  term  survival  of  the  CPSD  tiger  beetle  population  will  most  likely  occur  if  2.DD^ 

growth  rate  of  the  population  is  higher  than  the  mean  for  the  past  10  years  (which  ^ 

includes  7 years  of  negative  growth).  Management  options  to  increase  growth  rate  are 

limited  since  we  believe  that  climate,  especially  rainfall  and  soil  moisture,  is  the  main 

factor  affecting  growth  rate.  Although  carrying  capacity  increase  had  only  a limited 

effect,  increasing  the  amount  of  protected  area,  supplementing  moisture  in  localized 

patches,  or  reducing  vegetation  in  selected  swales  might  increase  available  habitat,  and 

should  be  considered  on  an  experimental  basis. 


PART  I.  POPULATION  AND  BIOLOGICAL  STUDIES 


INTRODUCTION 

The  Coral  Pink  Sand  Dunes  Tiger  Beetle,  Cicindela  albissima,  is  known  only 
from  the  Coral  Pink  Sand  Dunes  in  southern  Utah.  Studies  on  this  rare  insect  were  begun 
in  1991  when  a population  census  and  other  preliminary  studies  were  conducted  and 
suggested  its  probable  rarity.  Several  years  later,  a taxonomic  study  of  this  beetle  and  its 
relatives  in  the  Cicindela  martima  group  using  mitochondria  DNA  indicated  this  beetle 
was  sufficiently  distinct  from  the  other  subspecies  of  Cicindela  limbata  that  it  should  be 
considered  a separate  species.  Studies  on  its  biology,  especially  the  population  dynamics, 
have  been  carried  out  each  year  since  those  initial  studies.  These  studies  have 
documented  a highly  variable  population  size  for  C.  /.  albissima  since  1 992,  with 
numbers  ranging  from  a peak  of  2944  in  1996  to  a low  count  of  558  in  2005.  The  report 
presents  results  of  our  continuing  studies  with  this  rare  beetle  at  Coral  Pink  Sand  Dunes, 
with  support  from  a Utah  BLM  Challenge  Cost  Share  Program.  Specific  research 
objectives  for  the  2008  field  season  included  continued  monitoring  population  size 
throughout  the  dune  field,  studying  impact  of  rainfall  and  soil  moisture  on  population 
dynamics,  and  developing  a Population  Viability  Analysis  using  new  information 
obtained  since  1998. 


METHODS 


Adult  Surveys 

Throughout  the  years  of  this  study,  monitoring  adult  population  size  has  been  a 
primary  major  research  goal,  and  again  in  2008  this  was  continued  in  the  same  occupied 
habitat  of  the  dunes  using  the  same  method  as  in  previous  years.  Knowing  the  population 
size  in  each  year  is  important  for  documenting  the  trend  in  numbers  over  time  and  as  an 
indicator  of  the  long-term  viability  of  the  beetle  population  and  health  of  the  habitat  at 
CPSD.  The  adult  surveys  in  2008  were  done  from  May  13-24,  the  same  time  period  as  in 
previous  years.  The  method  we  used  to  determine  adult  size  in  2008  and  since  1999  was 
the  removal  method.  To  obtain  accurate  results,  the  survey  work  needs  to  be  carried  out 
during  periods  of  peak  adult  activity:  sunny  and  warm,  when  sand  surface  temps  are 
warm  enough  to  initiate  adult  activity  and  not  too  hot  to  be  unfavorable.  On  most  days. 
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these  optimum  conditions  for  conducting  the  removals  occurred  between  930  and  1230- 
1330  h and  again  from  about  1630-1830  h.  On  sunny  days  when  surface  temperatures 
increase  to  high  levels  (>45-50  C.),  adults  begin  digging  burrows  to  escape  the 
unfavorable  conditions.  They  typically  re-emerge  about  1630-1830  for  another  2 hour 
activity  period  before  digging  burrows  again  to  spend  the  night 

The  usual  procedure  for  conducting  the  removal  method  involves  three  workers 
systematically  moving  slowly  through  the  area  of  one  swale  and  collecting  all  of  the 
adults  that  are  observed.  Typically,  a small  number  of  the  adults  in  some  swales  are 
missed,  flew  off,  or  otherwise  not  captured.  We  have  estimated  that  this  loss  accounts  for 
only  about  5%  of  the  population.  After  completing  the  first  search  and  capture  in  a 
swale,  the  captured  adults  were  counted,  placed  in  individual  vials  and  retained  in  a 
cooler  with  ice.  We  then  repeated  the  process  a second,  third  and  sometimes  fourth  time 
until  we  captured  decreasing  numbers  with  each  pass.  Collected  adults  were  placed  in 
vials  and  stored  in  the  cooler  after  each  pass.  After  the  last  pass  all  of  the  adults  were 
released  throughout  the  swale  areas  where  they  were  collected.  We  then  moved  on  to 
another  swale  and  repeated  the  process  if  conditions  were  suitable.  For  many  of  the 
swales  that  had  few  adults  (<10-1 5)  we  made  only  one  pass  and  used  this  number  as  the 
estimate  for  that  swale.  The  removal  results  for  each  swale  were  analyzed  using  Program 
CAPTURE  to  produce  population  estimates  and  associated  variances  with  totals  from 
index  counts  added  in  to  produce  a total  estimate  for  the  population  and  the  95% 
confidence  intervals.  After  each  swale  was  sampled  we  walked  the  perimeter  of  the 
specific  area  sampled  with  a GPS  and  subsequently  produced  a GIS  map  of  the  whole 
survey  area.  From  these  maps  and  the  survey  data  we  computed  densities  of  adults  per 
swale  so  we  could  evaluate  the  habitat  quality  of  swales  for  supporting  beetles. 

Additional  surveys  using  the  visual  index  method  were  also  carried  out  in  several  swales 
in  area  AAA  and  AA  at  the  north  end  of  the  dune  field  and  in  area  D to  the  southwhere 
adults  have  been  foimd  in  previous  surveys. 

Rainfall 

In  this  report  we  include  monthly  rainfall  records  from  Kanab  for  the  months  that 
tiger  beetle  adults  and  larvae  are  active  (April  through  October)  for  the  years  1992 
through  2008.  We  use  the  Kanab  station  data  because  it  provides  the  most  complete 
record  and  because  the  rainfall  gauge  at  Coral  Pink  Sand  Dunes  State  Park  has  been 
periodically  non- functional.  Rainfall  and  the  coincidental  soil  moisture  are  believed  to  be 
the  most  important  factor  driving  the  population  dynamics  of  the  CPSD  tiger  beetle 
because  increased  rainfall/moisture  has  significant  positive  effects  on  many  life  history 
parameters,  including  oviposition,  recruitment,  survivor,  development  rate,  and  prey 
abundance. 

Larval  Surveys  and  Translocation 

We  did  not  conduct  larval  (or  adult)  translocations  in  either  May  or  September  of 
2008  because  of  the  small  population  size  of  adults  and  the  corresponding  small  numbers 
of  larvae. 
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Population  Viability  Analysis 


A major  part  of  the  2008  study  was  conducting  a population  viability  analysis  (PVA)  for 
the  CPSDTB.  This  analysis  was  previously  done  in  1999  based  on  1992-1998  population 
estimates  but  was  redone  because  much  additional  population  information  has  been 
obtained  and  our  population  estimates  since  that  time  have  been  more  accurate.  The 
analysis  used  the  same  basic  approach  as  the  previous  PVA  but  included  new  input 
parameters.  This  study  is  included  in  Part  2 of  this  report. 


RESULTS  AND  DISCUSSION 

The  total  adult  population  size  estimate  of  C.  albissima  for  2008  was  1072, 
with  a 95%  confidence  range  of  959  to  1185.  This  count  represents  a significant 
increase  fi-om  the  700  in  2007,  but  a decline  from  the  previous  even  year  cohort  estimate 
of  1 1 12  in  2006  (Fig.  1).  The  2008  estimate  also  indicates  a 6 year  period  of  counts  of 
1 1 12  or  less  since  the  highest  ever  count  of  2944  in  2002.  Moreover,  these  counts 
include  a much  higher  proportion  of  lowest  counts  for  any  time  since  counts  were 
initiated  in  1992.  Althou^  many  of  the  earlier  counts  (1992-1998)  were  determined  to 
be  overestimates,  the  accurate  estimates  from  1999-2001  (1804, 1270, 1459)  all  were 
higher  than  any  estimates  in  the  past  six  years.  We  have  attributed  these  population 
dynamics  and  the  recent  low  coxmts  to  the  significantly  reduced  rainfall  during  this  period 
(see  below). 

Fig.  1.  Population  estimates  of  CPSD  Tiger  Beetle  based  on  index  counts  (all  years) 
and  other  methods  (mark-recapture  method,  1992-1998  and  removal  method,  1999- 
2008).  Studies  have  indicated  values  for  1992-1998  overestimated  population  size 
while  those  from  1999-2008  are  reliable  estimates. 
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The  densities  and  distribution  pattern  of  adults  among  the  individual  swales  were 
generally  similar  to  most  of  the  previous  results  (Fig.  2,  Table  1).  Highest  adult  numbers 
were  in  swales  JK  with  280,  2-3  with  171,  H-WH  with  161,1  with  118,  publie  with  63 
and  QR/EQR  with  107.  As  has  been  the  norm  in  low  abundance  years,  numbers  of  adults 
in  swales  outside  of  the  primary  core  area  were  very  low.  Only  1 5 adults  were  found  in 
all  swales  north  of  swale  4 and  only  1 1 adults  were  in  the  D swales  which  are  south  of  the 
conservation  area.  Numbers  of  adults  and  larvae  also  generally  low  in  the  far  northern 
part  of  the  dune  field.  Adult  numbers  have  always  been  very  low,  typically  5-15  total  in 
this  area  and  only  2 were  counted  in  2008  (both  in  AAA-2).  Larval  numbers  in  these 
northern  swales  have  been  more  variable  and  higher  than  adult  numbers,  ranging  fi'om 
452  in  1999  to  0 in  2005  and  2006  (Fig.  3).  The  larval  count  in  2008  was  182,  160  of 
these  in  AAA-2;  2 first  instars,  52  seconds,  and  128  thirds.  This  relatively  high  and 
increasing  number  of  larvae  in  the  past  two  years  is  unusual  since  previously  higher 
numbers  in  the  northern  swales  have  been  associated  with  higher  rainfall  and  higher  adult 
numbers  in  the  main  population.  Possibly  there  may  have  been  more  dispersal  and/or 
higher  rainfall  in  this  area  in  2007,  which  would  account  for  the  higher  larval  numbers 
(especially  third  instars)  in  2008. 


Table  1.  Results  of 2008  Removal  Method  Estimates  and  Index  Counts  for  CPSD 
Tiger  Beetle  in  individual  swales.  N-hat  is  the  number  estimated  for  the  respective 
swales. 


2008  RESULTS 
THREE  PASS  SWALES 


Pass 


Swale 

1 

2 

3 

Total 

N-hat 

SE 

Var 

Public 

24 

15 

10 

49 

63 

11.207 

125.601 

JK 

178 

65 

24 

267 

280 

5.4263 

29.445 

E2&3 

13 

3 

1 

17 

17 

0.5312 

0.282 

E1&4 

12 

9 

9 

30 

60 

42.7066 

1823.854 

HWH 

94 

35 

20 

149 

161 

6.1359 

37.649 

1 

56 

28 

17 

101 

118 

9.2702 

85.937 

WY 

21 

11 

6 

38 

43 

4.6315 

21.451 

TOTAL 

398 

166 

87 

651 

742 

NA 

2124.219 

FOUR  PASS  SWALES 

Pass 

Swale 

1 

2 

3 

4 

Total 

N-hat 

SE 

Var 

2&3 

52 

54 

31 

16 

153 

171 

9.7979 

95.999 

QR 

22 

16 

16 

11 

65 

107 

33.1424 

1098.4 

0 

0 

TOTAL 

74 

70 

47 

27 

218 

278 

NA 

1194.4 

Total  beetles  handled  869 

Total  Estimate  1072 

Confidence  limits:  959  to1185 
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Swale 

N-hat 

SE 

Var 

Public 

63 

11.2072 

125.60133 

JK 

280 

5.4263 

29.444732 

E2&3 

17 

0.5312 

0.2821734 

E1&4 

60 

42.7066 

1823.8537 

HWH 

161 

6.1359 

37.649269 

1 

118 

9.2702 

85.936608 

WY 

43 

4.6315 

21.450792 

2&3 

171 

9.7979 

95.998844 

QR 

107 

33.1424 

1098.4187 

Total: 

1020 

3318.636 

Index  Count  Swales 

Swale 

Estimate 

LM 

1 

NA 

Index  Count 

N 

1 

NA 

Index  Count 

OP 

3 

NA 

Index  Count 

I 

2 

NA 

Index  Count 

DT 

3 

NA 

Index  Count 

E6&7 

2 

NA 

Index  Count 

4 

17 

NA 

Index  Count 

5 

7 

NA 

Index  Count 

9 

0 

NA 

Index  Count 

6&7 

5 

NA 

Index  Count 

8 

3 

NA 

Index  Count 

ST 

2 

NA 

Index  Count 

U 

0 

NA 

Index  Count 

V 

2 

NA 

Index  Count 

X 

2 

NA 

Index  Count 

Y 

1 

NA 

Index  Count 

WX 

1 

NA 

Index  Count 

52 

Fig.2.  Map  showing  location  of  individual  swales  at  CPSD.  Note  that  legend  is  a 
typical  density  pattern  and  not  specific  for  2008  results. 
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Fig.  3.  Numbers  of  larvae  counted  in  Conservation  Area  B (swale  AA-2  and  total  for 
all  swales)  at  the  north  end  of  the  CPSD  dune  field,  1996-2008. 


Adults  in  Protected  Versus  Unprotected  Swales 

Table  2 presents  the  numbers  of  adults  in  individual  protected  swales  (within  the 
conservation  area)  versus  numbers  in  the  unprotected  swales  outside  the  conservation 
area  (see  Fig.  2 for  locations).  The  relative  percents  in  these  two  areas  in  2008  were  very 
similar  to  the  values  of  most  recent  years.  In  2008,  80%  were  within  Conservation  Area 
A and  20%  were  in  the  unprotected  swales  outside  the  Conservation  Area.  Prior  to  the 
implementation  of  the  Conservation  Agreement  (1992  to  1998),  the  percent  in  the 
currently  protected  swales  were  70%  or  less.  This  pattern  of  a high  percent  within  the 
protected  swales  suggests  that  the  conservation  area  with  its  protected  swales  is  very 
important  in  supporting  the  core  of  the  population,  especially  during  years  when  the  total 
population  is  low.  The  Conservation  Area  not  only  protects  adults  and  larvae  from  OHV 
activity,  but  includes  the  highest  quality  swales  for  recruitment  and  adult  activity. 
Consequently,  this  protection  reserve  provides  the  best  conditions  for  protection  of  most 
of  the  core  population  and  can  thus  allow  for  higher  growth  rates  and  associated  dispersal 
to  other  areas  when  more  favorable  climatic  conditions  occur. 

Numbers  of  adults  have  been  variable  within  individual  swales  as  have  adults 
over  the  past  1 0 years,  but  some  swales  have  experienced  changes  independent  of  overall 
population  trends.  Swale  4 has  supported  many  few  adults  since  2001,  compared  to  years 
before  that  time.  The  only  notable  change  has  been  increasing  vegetation  and 
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coincidental  stabilization  of  portions  of  the  swale.  This  may  have  reduced 
oviposition/recruitment  areas  and  subsequently  adult  numbers.  Estimates  in  swales  2-3 
have  also  declined  hut  less  significantly,  perhaps  as  a result  of  a similar  increase  in 
vegetation  cover.  Among  the  unprotected  swales,  adult  numbers  have  exhibited  a pattern 
of  decline  in  swales  ST,  E-4,6,  E-6,7,  swales  N of  swale  4 (5-9),  OP  and  LMN. 
Vegetation  increase  cannot  explain  the  general  pattern  of  decline  in  these  swales,  hut 
perhaps  OHV  impact  or  other  factors  are  involved. 


Table  2.  Numbers  of  adults  in  swales  protected  within  Conservation  Area  A versus 
those  outside  this  area,  since  1998  when  protection  area  was  established. 


Swale  1998  1999 

2000 

2001 

2002 

2003 

2004 

2005 

2006 

2007 

2008 

CONSERVATION  AREA  SWALES 

4 

46 

98 

23 

62 

74 

20 

5 

9 

28 

4 

17 

1 

52 

142 

87 

116 

205 

48 

75 

72 

76 

45 

118 

pub 

76 

97 

275 

36 

65 

6 

155 

60 

63 

2,3 

155 

213 

244 

278 

349 

108 

111 

144 

148 

140 

171 

UK 

91 

266 

189 

191 

748 

114 

114 

91 

265 

143 

280 

H WH 

59 

369 

133 

180 

386 

90 

179 

111 

178 

114 

161 

W-WW 

31 

18 

3 

7 

15 

2 

1 

0 

0 

0 

0 

WX 

22 

66 

82 

42 
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Rainfall  Patterns  and  Effects  of  C.  albissima  Populations 

The  most  important  goal  of  our  studies  with  the  CPSD  Tiger  Beetle  is  to  obtain 
information  that  will  contribute  to  the  protection  and  recovery  of  this  rare  species. 
Understanding  what  factors  regulate  population  size  and  changing  dynamics  is  probably 
the  most  critical  aspect  of  this  goal.  As  a result  of  our  long  term  studies  with  this  beetle 
and  additional  experience  with  tiger  beetles,  we  have  become  convinced  that  rainfall  is 
the  primary  factor  controlling  population  size  and  the  changing  dynamics.  Included  in 
these  recent  studies  at  CPSD  are  field  experiments  that  have  documented  how 
supplemental  watering  patches  within  interdunal  swale  areas  attract  adults  with  a 
coincidental  increase  in  oviposition,  larval  recruitment,  larval  activity,  survival  and 
developmental  rate.  Evidence  has  also  demonstrated  how  increased  rainfall  and  soil 
moisture  increase  prey  abundance  which  contributes  to  increased  larval  survival  and  adult 
fecundity.  However,  our  attempts  to  directly  correlate  adult  numbers  with  rainfall 
amounts  over  the  past  1 5 years  have  not  been  as  conclusive.  A problem  in  developing 
this  evidence  is  the  two  year  life  cycle  of  the  tiger  beetle  and  the  varied  direct  and 
indirect  effects  of  rainfall/soil  moisture  on  many  life  history  parameters.  An  additional 
problem  may  be  our  use  of  rainfall  data  fi"om  the  Kanab  station  which  may  or  may  not  be 
closely  correlated  with  what  is  actually  occurring  at  the  dunes.  Kanab  is  the  closest 
station  available  for  use  and  should  provide  us  general  useful  trends  in  annual  rainfall 
over  the  years. 

There  is  some  evidence  that  our  rainfall  data  has  supported  the  positive  effects  of 
rainfall  on  the  adult  population  (Table  3,  Fig.  4).  For  example,  the  high  adult  numbers  in 

1996  and  2002  followed  several  prior  years  of  high  rainfall.  Also,  we  have  noted  reduced 
adult  numbers  in  2003-2005  following  low  rainfall  amounts  in  2001  to  2003.  Another 
complication  for  developing  correlations  is  that  the  actual  pattern  of  rainfall  within  a year 
may  be  more  important  than  the  total  rainfall.  For  example,  the  low  adult  numbers  in 
2004  and  2005  might  be  explained  by  the  low  rainfall  in  the  April-June  period  in  2002 
and  2003  compared  to  most  previous  years.  There  was  significant  rain  in  April  2003  but 
this  was  early  in  the  month  when  it  may  have  had  little  positive  effect  on  adult 
oviposition.  Also,  total  rainfall  (April-October)  for  2002  and  2003  were  the  lowest  of 
any  year  since  the  early  1990’s  and  with  low  spring  amounts  probably  combined  to 
explain  the  low  numbers  in  2003  and  2004.  In  contrast,  the  May  and  total  rainfall  from 

1997  to  2000  included  four  of  the  highest  rainfall  years  in  the  past  12  years  and  may  have 
explained  the  build-up  of  larvae  in  the  northern  dune  area.  Larvae  there  peaked  to  high 
numbers  in  1998  and  1999. 

The  rainfall  amounts,  both  total  and  April-June  in  2006,  2007  and  2008  have  been 
significantly  lower  than  any  other  three- year  period,  and  consequently  project  continuing 
low  adult  number  for  at  least  the  next  two  years.  The  lowest  ever  total  rainfall  for  Kanab 
was  recorded  in  2008  (2.82  inches)  and  perhaps  more  important  is  that  the  rainfall  for 
April-June  (0.22)  was  also  the  lowest  ever.  This  compares  with  the  previous  respective 
lows  of  5.15  in  1993  and  0.41  in  2002.  Rainfall  in  May,  probably  the  most  important 
month  for  recruitment,  was  only  .05  in  2006,  .04  in  2007,  and  .13  in  2008,  a much  lower 
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total  for  this  month  than  any  other  three  - year  period.  Our  research  plan  for  the  next  4 
or  more  years  is  to  put  in  place  a large  number  of  soil  moisture  probes  at  several  depths 
over  much  of  the  dune  field  to  track  weekly  and  seasonal  variations  in  soil  moisture. 


Fig.  4.  Rainfall  totals  for  April  through  June  and  April  through  October  at  Kanab, 
UT,  1992  to  2008.  April- June  is  the  critical  period  for  oviposition  and  larval 
recruitment 


—Annual  Total  —Apr-June  Total 
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Table  3.  Monthly  rainfall  at  Kanab  for  periods  when  adults  and  larvae  are  active, 
April-June  and  April-October,  1992-2008. 
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II.  POPULATION  VIABILITY  ANALYSIS  OF  C ALBISSIMA  AND 
IMPLICATIONS  FOR  MANAGEMENT 


Introduction 

Population  viability  analyses  are  an  increasingly  important  tool  for  studies  of  the 
conservation  of  rare  species.  PVA  models  have  been  used  to  assist  in  identifying  the 
threats  faced  by  a species  and  evaluating  the  likelihood  that  the  species  will  persist  for  a 
given  time  into  the  future.  A PVA  may  be  used  to  address  three  aspects  of  threatened 
species  management:  1 . It  may  reveal  that  a population’s  long  term  viability  is  insensitive 
to  certain  population  parameters  and  research  may  be  directed  at  parameters  that  do  have 
an  impact  on  potential  extinction;  2.  It  may  be  used  to  assess  the  vulnerability  of  a 
species  to  extinction,  and  consequently,  results  of  the  analysis  may  be  used  to  set  policies 
and  priorities  by  conservation  groups;  and  3.  It  may  be  used  to  predict  potential  effects  of 
reintroduction,  captive  rearing,  habitat  manipulation,  and  other  management  strategies. 

This  part  of  the  2008  Annual  Report  presents  a second  population  viability 
analysis  (PVA)  for  the  Coral  Pink  Sand  Dunes  Tiger  Beetle.  This  updated  version  of  the 
PVA  includes  the  same  basic  model  structure  and  the  same  parameters  as  the  1999 
version,  but  is  based  on  new  parameter  values  which  have  been  obtained  by  the 
population  studies  in  the  past  10  years,  1999-2008.  The  first  PVA  was  based  on  data  on 
the  population  from  1992  to  1998.  Since  the  first  PVA,  we  have  used  population 
estimates  that  are  more  accurate  and  thus  more  reliable  for  determining  inputs  for  this 
revised  PVA.  These  recent  studies  have  determined  that  the  population  estimates  used  in 
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that  previous  PVA  overestimated  adult  numbers  in  many  of  those  years  and  were  thus 
less  reliable.  Also,  in  part  because  of  these  new  estimates  and  also  because  the  true 
population  size  has  been  lower  in  these  past  10  years  than  in  the  previous  7 years,  a new 
model  is  very  important  for  re-evaluation  of  extinction  probabilities  which  will  thus 
provide  valuable  input  for  the  revision  to  the  Conservation  Agreement. 

METHODS 


Background 

PVA’s  are  usually  developed  from  age-  or  stage-based  population  models 
requiring  detailed  and  accurate  input  data,  such  as  transition  and  survival  probabilities 
between  life  stages  and  stage  specific  fecundity  rates.  Because  much  of  this  information 
is  not  available  for  the  CPSD  tiger  beetle,  an  alternative  method  based  on  empirical 
estimates  of  adult  population  size  obtained  each  year  from  1999-2008.  In  many  regards 
this  method  is  quite  appropriate  to  use  and  may  provide  a more  accurate  assessment  than 
stage-based  models  because  our  population  estimates  were  determined  to  be  very 
accurate  and  the  cumulative  errors  that  may  accrue  from  the  multiple  parameters  in  these 
other  models  are  absent. 

Modeling  Approach 

The  CPSD  tiger  beetle  population  was  assumed  here  to  consist  of  one  panmictic 
population,  based  on  field  observations  over  the  years  of  study,  1992-2008.  A small 
satellite  population  is  found  in  the  northern  part  of  the  dune  field,  but  is  considered  to  be 
a “sink”  that  contains  excess  individuals  that  have  dispersed  from  the  main  southern 
population  concentrated  within  Conservation  Area  B during  periods  of  high  density. 
Consequently,  the  model  did  not  use  a metapopulation  approach  which  is  common  to 
many  other  tiger  beetle  species  which  have  multiple  populations. 

In  selecting  a population  model,  we  analyzed  the  observed  population  dynamics 
of  the  tiger  beetle  population  in  its  primary  habitat  and  calculated  the  key  population 
parameters  to  be  used  in  the  model.  As  indicated  above,  unlike  the  previous  PVA,  we 
based  our  input  parameters  on  aimual  population  estimates  from  1999  through  2008. 
These  estimates  were  obtained  using  a removal  method  and  program  Capture  which 
consistently  provided  low  confidence  intervals  and  thus  accurate  estimates  of  the 
populations  each  of  these  years.  Mean  population  size  was  calculated  to  be  1227 
individuals,  with  total  population  estimate  being  as  low  as  558  and  as  high  as  2944  over 
the  10  years  (Table  1).  These  values  were  analyzed  to  determine  if  the  population 
exhibited  density  dependent  fluctuations.  Because  the  beetles  exhibit  a two-year  life 
cycle,  density  dependence  would  occur  on  a two  year  interval.  If  density  dependent 
effects  influenced  population  dynamics,  the  slope  of  a regression  of  population 
abundance  in  year  t (Nt  log-transformed)  on  the  rate  of  growth  between  year  t and  year  t 
+■7  = lambda,  log  transformed)  should  be  negative.  The  slope  was  not 

significantly  different  than  zero,  indicating  a lack  of  density  dependence.  However,  the 
statistical  power  of  the  test  was  low  due  to  the  limited  total  number  of  years  in  the 
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analysis  (8  cohorts),  and  it  is  possible  that  density  dependent  effects  operate  in  this 
population.  Morever,  it  is  apparent  from  experience  with  this  species  that  limitations 
related  to  the  amount  and  quality  of  habitat  (due  to  climatic  induced  variations  in  food 
availability,  OHV  disturbance  and  other  factors)  do  exert  an  influence  on  the  population. 
Because  of  this,  a logistic  growth  model  with  a ceiling  carrying  capacity  was  selected  for 
the  modeling. 

Selected  parameters 

The  calculated  parameters  used  as  the  basic  input  parameters  for  the  model 
based  on  the  1999-2008  population  data  and  additional  parameters  included  in  various 
model  runs  are  given  in  Table  2.  As  indicated  above,  the  mean  population  size  for  the  10 
year  period  was  1227,  and  thus  rounded  to  1200  for  the  starting  population  size  for  all 
runs.  The  calculated  mean  growth  rate  was  1.058  so  this  was  used  as  the  basic  value. 
Although  growth  rates  of  less  than  1.0  were  recorded  for  7 of  the  10  years,  indicating  a 
decline  from  one  year  to  two  years  later  (two  year  life  cycle),  we  used  only  positive 
values  of  growth  rate  in  the  model  runs  since  negative  values  would  result  in  relatively 
rapid  extinction.  This  is  unrealistic  for  the  history  of  this  population  which  has 
obviously  existed  for  a very  long  period  of  time  in  these  dunes  with  erratic  population 
fluctuations,  so  the  growth  rate  has  had  to  be  on  average  over  1.0.  Variation  in  growth 
rates  is  characteristic  of  insect  populations  which  exhibit  explosive  growth  during  periods 
of  favorable  conditions  and  dramatic  declines  when  conditions  are  unfavorable.  In  the 
various  model  runs  we  used  the  maximum  growth  rate  that  we  observed  (2.3)  and  also  an 
intermediate  value  of  1.35  which  we  believe  to  be  a value  that  is  realistic  for  this 
population. 

Table  1.  Population  estimates  and  calculated  cohort  growth  rates  of  the  CPSD  Tiger 
Beetle  Population,  1999-2008. 


Population 

Growth 

Year 

Estimate 

Rate 

1999 

1804 

0.81 

2000 

1270 

2.32 

2001 

1459 

0.41 

2002 

2944 

0.3 

2003 

595 

0.94 

2004 

757 

1.47 

2005 

558 

1.25 

2006 

1112 

0.96 

2007 

700 

2008 

1072 

mean 

1227 
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One  of  the  most  important  aspeets  in  PVA  is  to  aeeount  for  the  effeets  of 
stoehastie  events  on  population  growth  rate.  This  is  espeeially  true  for  invertebrate 
populations  wherein  temporal  variation  in  abundanee  is  probably  due  more  to  stoehastie 
events  in  the  environment  than  to  density  dependent  effects.  Environmental  stochasticity 
was  incorporated  into  the  model  by  making  cohort  population  growth  rate  a random 
variable.  For  each  specified  lambda  (growth  rate)  the  cohort  growth  rate  for  a particular 
time  step  in  the  model  was  selected  from  a normal  distribution  with  mean  growth  rate  and 
standard  deviation  equal  to  the  value  specified  by  the  user.  Larger  specified  values  of 
the  standard  deviation  in  the  model  mean  greater  potential  for  any  given  year  to  be 
extremely  poor  or  extremely  good  in  terms  of  population  growth  rate.  Appropriate  values 
for  the  standard  deviation  of  the  growth  rate  were  based  on  observed  variation  in  cohort 
growth  rates  between  1999-2008.  The  standard  deviation  in  cohort  growth  rate  was 
calculated  on  an  annual  basis  (i.e.,  as  if  odd  and  even  year  cohorts  did  not  exist).  This 
value,  .646,  was  used  as  the  base  value  in  the  model,  but  to  model  a broader  range  of 
potential  stochasticity  over  time,  we  selected  a lower  value  of  .33  and  a higher  value  of 
.95  as  additional  parameter  values  in  various  model  runs. 

The  effects  of  carrying  capacity  (K)  on  model  predictions  were  examined  by 
varying  K in  the  model  (K  = 800,  1200,  2000,  3500,  5000),  while  holding  other 
parameters  constant.  The  base  value  for  carrying  capacity  in  the  model  was  set  at  3500,  a 
value  that  was  slightly  higher  than  the  maximum  population  size  we  have  ever  recorded, 
2944.  We  increased  this  to  3500  to  account  for  some  additional  adults  that  would  likely 
be  missed  in  any  estimate,  as  we  believe  to  be  the  case. 


Table  2.  Parameter  values  used  in  this  PVA  and  those  used  for  the  previous  PVA. 


Parameters  For  Model 

Base  Values 
from  1999-2006  Data 

Alternative  Values  Used 

Base  Values 
1992-1998  Data 

Mean  Growth  Rate 

1.058 

1.35 

23 

.86/1.5* 

Std  Dev. 

0.646 

0.33 

0.95 

0.95 

Carrying  Capacity 

3500 

800 

1200 

2000 

5000 

3000 

Starting  Population 

1200 

2025 

Ramas  Models  Used 

Logistic  with  ceiling 

Saamble 

Contest 

Lx)g Stic  with  ceiling 

*1.5  was  selected  as  the  most  likely  growth  rate  in  the  previous  PVA.  The  calculated  mean  of  .86 
is  a negative  growth  rate  which  would  have  resulted  in  rapid  extinction. 


Predicting  Probability  of  Extinction  Under  Various  Model  Scenarios. 

Program  ECOLAB  (Ramas  Software,  Inc.)  was  used  to  model  population 
d5mamics  and  extinction  probability.  The  model  ran  for  50  time  steps  (equal  to  50  cohorts 
and  thus  100  years).  At  the  beginning  of  each  model  run,  the  starting  population  size  was 
specified  as  1200,  the  average  population  size  for  the  past  10  years.  Population 
abundance  for  the  next  time  step  was  predicted  based  on  the  model’s  predicted  lambda. 


16 


generated  from  the  normal  distribution  with  specified  standard  deviation.  This  predicted 
population  was  then  used  in  the  next  time  step.  The  process  was  repeated  for  50  time 
steps  producing  a time  series  of  abundance  for  a simulated  100  year  period.  When  a 
simulation  was  complete,  original  parameter  values  were  reset  to  the  specified  initial 
conditions,  and  the  model  rerun.  Predictions  from  each  individual  run,  even  if  started 
from  the  same  initial  conditions,  produced  different  abxmdance  time  series  because  of  the 
environmental  stochasticity  incorporated  into  the  model.  In  order  to  predict  extinction 
probability  under  a specified  set  of  initial  conditions,  1000  replicate  model  runs  were 
made.  The  percentage  of  runs  that  produced  an  extinction  (population  size  reaching  0 at 
some  point  during  the  100  year  simulation)  was  the  predicted  probability  of  extinction  for 
CPSD  tiger  beetles  over  the  next  100  years.  By  varying  initial  conditions  as  indicated  in 
Table  2,  the  effect  of  different  model  assumptions  on  predicted  probability  of  extinction 
was  evaluated. 

It  should  also  be  pointed  out  that  extinction  of  a population  might  also  occur  even 
before  the  population  reached  zero,  as  a result  of  Allee  effects.  Alice  effects  describe  a 
situation  where  a population  growth  rate  drops  steeply,  often  to  extinction,  once  a 
critically  small  population  size  is  reached.  For  example,  with  too  few  individuals,  finding 
mates  may  become  difficult,  leading  to  rapid  decline  in  abundance.  There  is  also 
evidence  in  many  species  that  low  population  size  of  as  little  as  500  to  50  produces 
genetic  decline  due  to  loss  of  genetic  diversity,  inbreeding  or  related  factors.  Small 
populations  are  also  more  susceptible  to  extinction  from  catastrophic  weather  events  or 
other  natural  or  man-induced  disturbances.  Consequently,  our  results  in  all  model  runs 
show  extinction  risks  to  0,  50,  100,  200  and  500  individuals. 


RESULTS  AND  DISCUSSION 

The  results  of  the  series  of  model  runs  using  the  range  of  parameters  produced 
extinction  rates  from  zero  (no  predicted  chance  of  population  extinction)  to  1.0  (a  100% 
chance  of  extinction).  Not  surprising  was  that  varying  model  parameters  had  dramatic 
effects  on  predicted  extinction  probability  (see  Table  3).  This  was  determined  by  a 
sensitivity  analysis  in  which  one  parameter  was  varied  while  holding  the  others  constant. 
Several  parameters,  carrying  capacity  and  starting  population  size  had  relatively  little 
effect  on  extinction  probability.  For  example,  using  the  base  values  that  we  calculated 
for  growth  rate  and  standard  deviation,  probability  of  extinction  was  determined  to  be  .32 
(=32%  extinction  probability)  with  a 3500  carrying  capacity  (K),  .3 1 with  a K of  2000, 

.35  if  K was  1200,  and  .37  if  K was  800  (Table  3).  Starting  population  size  had  even  less 
affect  although  we  do  not  present  results  here.  Using  two  other  models  (logistic  scramble 
and  logistic  contest)  also  did  not  produce  appreciably  different  results  (Table  3). 

The  two  parameters  which  had  very  significant  effects  on  extinction  probabilities 
were  growth  rate  and  standard  deviation  (stochasticity)  (Table  3,  Fig.  1).  For  example, 
extinction  probability  is  zero  with  growth  rates  of  1.35  and  2.3  (when  stochasticity  is 
.65)  but  .32  when  growth  rate  is  1 .06.  Extinction  risk  increases  from  zero  to  over  .90 
when  stochasticity  increases  from  .33  to  .95.  The  results  of  all  combinations  included  in 
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the  PVA  model  runs  are  given  in  Table  3.  In  general  these  runs  show  that  the  CPSD  tiger 
beetle  population  is  at  risk  of  extinetion  if  growth  rate  remains  at  1 .058,  the  average  for 
the  past  8 eohorts),  even  when  stoehastieity  is  lower  (.33)  than  we  have  reeorded.  By 
eontrast,  growth  rates  of  1 .35  greatly  reduee  extinetion  risk,  unless  stoehastieity  is  .95. 

A growth  rate  of  2.3,  the  highest  reeorded  for  any  eohort,  makes  extinction  unlikely, 
even  when  stoehastieity  is  .95  (Table  3).  The  scenarios  which  we  believe  provides  the 
“best  guess”  of  a range  of  extinction  risks  are  the  one  with  the  calculate  parameters  for 
1999-2008  (1200  starting  population,  3500  carrying  capacity,  .646  standard  deviation  or 
stochastacity  measure,  and  1.058  growth  rate)  and  the  other  scenario  with  the  same 
parameters  except  a growth  rate  of  1 .35.  With  a 1 .058  growth  rate,  the  probability  of  the 
population  going  to  0 was  .32,  a high  likelihood  of  extinction,  but  a 1 .35  growth  rate 
produced  an  extinction  probability  of  only  .001 . Additional  output  from  the  most  likely 
parameter  values  in  model  runs  that  show  population  numbers  over  time  and  the  full 
range  of  extinction  probabilities  are  given  in  Fig.  2. 


Fig.  1.  Extinction  risks  resulting  from  different  parameters.  Probabilities  that  the 
population  will  reach  0,  50  or  100. 

A.  Extinction  probabilities  based  on  growth  rate  variations  (.65, 1.3,  2.3) 
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B.  Extinction  probabilities  based  on  stochasticity  (standard  deviation)  variations. 
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Table  3.  Results  of  model  runs  using  a range  of parameter  manipulations.  Albl,  etc 
are  designations  for  the  respective  parameters  used.  A starting  population  size  of  1200 
used  in  all  runs.  Note  thatAlbl4  is  comparable  to  the  “best  guess”  input  parameters 
for  the  previous  PVA. 


Albl 

Alb2 

Alb3 

Alb4 

Alb5 

Alb14 

Alb6 

Alb7 

Alb8 

Albl  7 

Albl  8 

Mean  Growth  Rate 

1.058 

1.35 

2.3 

1.058 

1.058 

1.35 

1.058 

1.058 

1.058 

1.058 

1.35 

Std  Dev. 

0.646 

0.646 

0.646 

0.33 

0.95 

0.95 

0.646 

0.646 

0.646 

0.646 

0.646 

Canving  Capacity 

3500 

3500 

3500 

3500 

3500 

3500 

800 

1200 

2000 

5000 

5000 

Prob.  Popul  will  reach: 

0 

0.32 

0.001 

0.0001 

0.002 

0.95 

0.007 

0.37 

0.35 

0.31 

0.27 

0.001 

50 

0.87 

0.003 

0 

0.09 

0.99 

0.22 

0.95 

0.94 

0.91 

0.84 

0.002 

100 

0.87 

0.002 

0 

0.18 

0.99 

0.36 

0.95 

0.94 

0.91 

0.91 

0.009 

200 

0.97 

.04 

0 

0.33 

0.99 

1 

1 

0.99 

0.96 

0.96 

0.02 

500 

0.99 

0.23 

0 

0.63 

1 

1 

1 

1 

1 

0.99 

0.16 

Albll 

Alb  15 

Alb16 

Albl  7 

Alb9 

AlblO 

Alb12 

Albl  3 

Albl  cor 

Alt2cxn 

Albl  sc 

Alb2sc 

Alb2sc 

Mean  Growth  Rate 

1.35 

23 

23 

1.058 

1.35 

1.35 

23 

1.058 

1.058 

1.35 

1.058 

1.35 

135 

Std  Dev. 

0.33 

0.96 

a33 

6646 

0.646 

633 

0.646 

0.95 

0.646 

0.646 

0.646 

6646 

0.646 

Canving  Capacity 

3500 

3500 

3500 

5000 

1200 

1200 

1200 

1200 

3500 

3500 

3500 

3500 

3500 

Prob.  Popit  mil  reach: 

C 

0 

0 

0 

0.27 

0.001 

0 

0 

0.88 

0.32 

0.001 

629 

6001 

0.001 

50 

0 

0 

0 

0.84 

0.012 

0 

0 

0.99 

69 

603 

-.99 

0.05 

0.05 

100 

0 

0 

0 

0.91 

0.055 

0 

0 

1 

0.94 

612 

694 

0.13 

0.13 

200 

0 

0 

0 

0.96 

0.213 

0 

0 

1 

0.98 

636 

0.97 

0.35 

035 

500 

aool 

0.001 

0 

0.99 

693 

0.002 

0.001 

1 

0.99 

688 

699 

0.87 
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Abundance 


Fig.  2.  Additional  graphical  representation  of  selected  model  runs. 

A.  Trajectory  graph  of population  size  over  time  for  Alb  1(1.058,  .646,  3500) 
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B.  Probabilities  population  will  reach  a particular  population  size  for  Alb  1 (1.058, 
.646,3500). 
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C.  Trajectory  graph  of population  size  over  time  for  Alb  2(1.35, . 646,  3500) 


Trajectory  summary 


D.  Probabilities  population  will  reach  a particular  population  size  for  Alb  2 (1.35, 
.646,  3500). 


Extinction/Decline 
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Abundance 


E.  Trajectory  graph  of population  size  over  time  for  Alb  14(1.35,  .95,  3500) 


Trajectory  summary 


F.  Probabilities  population  will  reach  a particular  population  size  forh  Alb  2 (1.35, 
.646,  3500). 
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SUMMARY  AND  M ANAGEMENT  CONSIDERATIONS 


The  results  of  this  updated  PVA  indieates  the  CPSD  tiger  beetle  is  probably  at 
even  greater  risk  of  extinction  that  projected  in  the  1999  PVA.  The  increased  risk  is 
essentially  a result  of  the  lower  population  sizes  recorded  in  the  past  1 0 years  and  the 
associated  low  annual  growth  rates  during  this  period  that  drove  the  model  results.  The 
low  growth  rates  resulting  from  low  recruitment  and  survival  during  this  1 0-year  period 
are  likely  a result  of  low  rainfall,  especially  in  the  past  5-6  years.  Unfortunately, 
management  options  for  increasing  population  size  are  limited  because  of  the  dominating 
effect  of  growth  rates  on  population  size.  Since  rainfall/soil  moisture  seems  to  be  so 
important,  supplementing  moisture  may  offer  an  option  for  increasing  ovipostion  and 
recruitment,  at  least  on  a limited  or  localized  basis.  Our  previous  experiments 
demonstrated  positive  effects  of  moisture  supplementation  in  localized  areas.  It  is 
reasonable  to  conclude  that  larger  scale  moisture  supplementation  during  the  critical 
period  of  adult  oviposition  (May)  in  selected  swales  could  have  measurable  effects,  and 
will  be  carried  out  in  2009.  A second  management  option  that  might  have  positive 
effects  on  growth  rate  by  increasing  oviposition  and  recruitment  is  to  reduce  vegetation 
cover  in  selected  swales  where  it  may  have  reduced  oviposition  sites  in  recent  years. 

This  option  might  also  be  implemented  on  a limited  scale  in  2009  by  manual  removal  of 
vegetation  (by  clipping)  is  several  swales.  A third  option  that  should  be  considered  is  an 
expansion  of  Conservation  A to  include  several  important  swales  that  are  believed  to 
have  suitable  habitat,  but  are  being  impacted  by  heavy  OHV  use  which  prevents 
successful  colonization  and  recruitment.  One  additional  management  option  to  increase 
growth  rate  and/or  reduce  extinction  risk  is  to  conduct  captive  rearing  trials  to  produce 
first  instar  larvae  for  translocation  into  Conservation  B in  an  attempt  to  establish  an 
additional  population  there. 
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